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The molecular behaviors of interfacial water molecules at the solid/liquid interface are of a fundamental 
significance in a diverse set of technical and scientific contexts, thus have drawn extensive attentions. On 
certain surfaces, the water monolayer may exhibit an ordered feature, which may result in the novel wetting 
phenomenon. In this article, based on the molecular dynamics simulations, we make a detailed structure analysis 
of the ordered water monolayer on ionic model surface with graphene-like hexagonal lattices under various 
charges and unit cell sizes. We carefully analyze the water density profiles and potential of mean force, which 
are the origin of the special hexagonal ordered water structures near the solid surface. The number of hydrogen 
bonds of the ordered water monolayer near the solid surface is carefully investigated. 
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I. INTRODUCTION 


The complex behaviors of interfacial water [1-9], which 
are of great importance in research fields of protein stability 
and folding [10], molecular self-assembly [11], manipulat- 
ing biomolecules [12], rearrangement of immunodeficiency 
virus [13] etc., have drawn extensive attentions [5—7, 14], 
since the molecular structure and dynamics of the interfacial 
water molecules are usually different from the bulk proper- 
ties [15]. Interfacial water molecules play an important role 
in biophysical process. For example, water effectively cat- 
alyzes chiral interconversion of thalidomide [16], and dewet- 
ting transition promotes the amyloid fibrils formation [17]. 
Owing to the interaction between the interfacial water and 
the hydrophilic solid substrate, the diffusion of interfacial wa- 
ter [18] is slower, and the lifetime of hydrogen bonds [19] is 
longer, than that of the bulk water, as having been confirmed 
by experiments [20-22]. Recently, ordered structure of the 
interfacial water confined [23] at one or two dimensions has 
been studied extensively by both experimental and theoret- 
ical methods. In 2009, we reported a liquid water droplet 
on a water monolayer, termed as “ordered water monolayer 
does not completely wet water” on a model surface at room 
temperature [24]. Later, similar phenomena were observed 
by several experiments on sapphire c-plane electrolyte sur- 
face [25] and on self-assemble monolayer (SAM) surfaces 
with the —COOH terminal [26, 27]. In addition, theoreti- 
cal simulations found similar phenomenon on hydroxylated 
metal oxide surfaces of Al,O3 and SiO, [4], Talc [28] and 
Pt(100) metal surfaces [29]. We also explored the effect of 
morphology [30] and the critical length of the charge dipoles 
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of the solid surface [31] on the structures of interfacial water 
and the surface wetting behaviors. 

In this article, based on molecular dynamics simulations, 
we investigate the structure and hydrogen bonds to show 
detail information of the ordered water monolayer on ionic 
model surface having graphene-like hexagonal lattices with 
various charges and unit cell sizes. The article is organized as 
follows. The ordered structure of water monolayer near the 
surface is described in Sec. III. A. In Sec. II. B, the water 
density and the potential of mean force (PMF) [28] are stud- 
ied. In Sec. III. C, the number of hydrogen bonds is calculated 
to show the stable formation of hydrogen bonds network in 
the ordered water monolayer. Finally, a short conclusion is 
presented in the last section. 


Il. SIMULATION DETAILS 


We configured a hexagonal solid lattice with 1664 solid 
atoms and the neighbor bond length / was described in Fig. 1, 
the same as our previous studies [24, 32]. The initial sys- 
tems for the molecular dynamics simulations contained a wa- 
ter layer of about 4.0 nm thick on the ionic model surface, 
where positive and negative charges were located diagonally 
in neighboring hexagon, and it was found that the charge 
had great influence on the flux of water molecules in nan- 
otube [12, 32]. All the simulations were performed at T = 
300 K (NVT ensemble), with Gromacs 4.5.4 [34] by using a 
time step of 1.0 fs. The Lennard-Jones parameters of the solid 
atoms were €,, = 0.105kcal/mol and css = 33.343 Å, and 
SPC/E water model [35] was used. The particle-mesh Ewald 
method [13] with a real space cutoff of 1 nm was adopted for 
the long-range electrostatic interactions and a 10 A cutoff was 
used for the van der Waals interactions. The periodic bound- 
ary conditions were applied in three directions. The simula- 
tion time for every system was 4 ns and the last 2 ns data was 
collected for analysis. 

Two series of simulations were performed to investigate 
the ordered water monolayer formation on a hexagonal 
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Fig. 1. (Color online) Hexagonal solid surface with charged pairs (a) and snapshot of the 4.0 nm thick water layer on the surface (b). The red 
and green spheres represent atoms with positive and negative charges, respectively, while neutral solid atoms are shown in white. 


polarity solid surface. In the first series of simulations, 
the charge q of the solid atoms increased from 0.6e 
to 1.0e with 0.le interval and there were 5252 water 
molecules in the simulation boxes with the volume of 
6.395nm x 6.816nm x 20.110nm. The value of the neigh- 
boring bond length of solid atoms was kept as the constant 
of 1 = 0.142nm. In the second series of simulations, the 
bond length / was set at 0.120nm, 0.130nm, 0.142 nm, 
0.150 nm and 0.160 nm, with q = 0.8e, and the water layer 
thickness was kept at about 4.0 nm, with the water molecules 
of 3525, 4314, 5252, 5721 and 6564, in the simulation 
boxes of 5.404nm x 5.760nm x 20.110nm, 5.854 nm x 
6.240 nm x 20.110 nm, 6.395 nm x 6.816nm x 20.110 nm, 
6.755nm x 7.20nm x 20.110nm and 7.205 nm x 
7.680 nm x 20.110 nm, respectively. 


HI. RESULTS AND DISCUSSIONS 
A. Structure analysis of the water monolayer 


To study the structure of water molecules in the water 
monolayer on the solid surface, two angle parameters 0 and 
ọ are introduced as illustrated in Fig. 2(c) and 2(d), where 
0 is defined as the angle between a water molecule dipole 
and z axis, and y is the angle formed between the projection 
onto x-y plane of a water dipole and a crystallographic di- 
rection [30]. Here, the definition of the water monolayer is 
the water molecules in the first layer next to the solid surface 
with an average thickness of 0.4nm, the same as our previous 
work [24], which is also consistent with the existence of an 
experimentally observable monolayer [36]. The second layer 
is defined as the water molecules with an average thickness 
of 0.4nm above the water monolayer. 

As shown in Fig. 2(e), two peaks of angle 0 confirm the 
two states, namely, state | and state 2 as depicted in Figs. 2(a) 
and 2(b). The left peaks at 0 ~ 60° represent state 1 with oxy- 
gen atoms attracted by the positive charged atoms, while the 


right peaks at 9 ~ 120° represent state 2 with —OH bonds 
pointing towards the negative charged atoms. Fig. 2(f) is 
the normalized probability distributions of angle y with three 
peaks at y = 0°, 120° and 240°, which demonstrate that the 
water molecules in the monolayer can form a 2D hexagonal 
configuration (Fig. 2(d)), the same as our previous work [24]. 
As q increases from 0.6e to 1.0e, all the peaks in Figs. 2(e) 
and 2(f) become higher and the water molecules in the mono- 
layer become ordered due to the larger binding of the surface 
charged atoms. However, the peaks are quite different as the 
bond length / increases. At l = 0.142 nm and q = 0.8e, the 
peaks are the highest (Figs. 2(g) and 2(h)), hence the most 
ordered water molecules in the water monolayer. As the l de- 
parts from 0.142 nm, the ordered hexagonal water monolayer 
gradually disappears. These results show that the ordered wa- 
ter structure greatly depends on the surface charge and suit- 
able cell size. 


B. Water density distribution profiles and PMF curves 


Figures. 3(a) and 3(b) show the water density as a function 
of z at different q and l. The reference z = 0 corresponds 
to the solid surface. Two peaks can be seen for all curves lo- 
cating at z = 0.3nm and 0.6nm. Due to the strong binding 
of charges on the surface, we can observe a quite high den- 
sity peak near the solid surface forming the monolayer. The 
density increases with the charge, reaching the largest at q = 
0.9e and 1.0e. With increasing cell size, the density increases 
first until 7 = 0.142 nm, where it begins to decrease, indicat- 
ing the formation and break-down of the ordered structure, 
respectively. 

The density relates to potential of mean force (PMF), F (z), 
by the expression [28], 


F(z) = —kpT In(p(z)/pw), (1) 


where, kg is the Boltzmann constant and p,, = 33nm~? is 


the number density of bulk water. F(z) is the potential of 
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Fig. 2. (Color online) (a) State 1 with water molecule adsorbed by positive binding charge and three negative neighbor charges marked with 
A, A’, B. (b) State 2 with water molecule adsorbed by negative binding charge and three positive neighbor charges marked with C, D, D’. (c) 
Schematic of angle 0 defined as the angle between a water molecule dipole and z axis. (d) Schematic of angle y defined as the angle formed 
between the projection onto x-y plane of a water dipole and a crystallographic. (e) Probability distribution of 0 in the monolayer vs. q. (f) 
Probability profile for angle y in the monolayer vs. q. (g) Probability profiles of 0 in the monolayer vs. l. (h) Probability distribution of y in 


the monolayer vs. l. 


mean force for bringing a water molecule from the bulk to a 
distance z from the solid surface. Figs. 3(c) and 3(d) show the 
PMF curves and for every curve there are two valleys at z = 
0.3nm and 0.6nm. The two valleys account for the adsorp- 
tion of the solid surface. The minimum PMF at z = 0.3nm 
is about —0.9 kcal/mol at q > 0.8e and! = 0.142nm. The 


PMF reveals the adsorption interaction of the solid surfaces 
at the valleys. The adsorption increases with the charge, dis- 
playing a wide range of binding strength to attract the water 
molecules and form the ordered monolayer. This is different 
from the bulk water. Suitable cell size is quite important for 
adsorption interaction of the solid surface and formation of 
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Fig. 3. (Color online) (a) Density profile of water molecules away from the surface vs. q, divided by the number density of bulk water, py = 
33nm~*. (b) The density profile vs. l, pz /pw. (c) Potential of mean force F (z) vs. q. (d) F(z) vs. L. 
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Fig. 4. (Color online) Average number of hydrogen bonds of a water molecule to other water molecules in the same layer (¢), to water 
molecules in the second layer (e), and their sum (W) as function of q (a) and / (b). 


the monolayer. The PMF results indicate that the distribution 
of water molecules and formation of the ordered water mono- 
layer are affected by the charge and cell size. 


C. Hydrogen bonds in the water monolayer 


The ordered water monolayer affects the formation of hy- 
drogen bonds of the water molecules in the interface. We 
calculated the average hydrogen bonds of a water molecule to 
its neighboring water molecules in the same monolayer (“in 
the monolayer” H bonds), and to water molecules in the sec- 
ond layer (“to the second layer” H bonds), as shown in Fig. 4. 


The criteria characterizing existence of hydrogen bond be- 
tween two water molecules is the geometric definition that 
their O—O distance is less than 3.5 Å and simultaneously the 
angle H—O - - - O is less than 30° [37]. 


In Fig. 4(a), the number of hydrogen bonds within the 
monolayer increases and the number of the hydrogen bonds 
between the monolayer and the second layer decreases as the 
increase of charge. Their sum remains at ~ 2.9 when q > 
0.8e, which approaches 3, the maximum number of hydro- 
gen bonds that any water molecule can form in the mono- 
layer [24]. The interaction energy between the monolayer 
and the charged surface is stronger when the charge increases 
as we calculate in Sec.II(B). The water molecules bound in 
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the monolayer make it easy to form hydrogen bond with the 
water molecule in the same layer. There is competition for 
formation of hydrogen bonds between the “in the monolayer” 
H bonds and “to the second layer” H bonds. The increase 
of former leads to the decrease of latter for weaker interac- 
tion between the water molecules in the monolayer and wa- 
ter molecules in the second layer. In Fig. 4(b), when / = 
0.142nm and q = 0.8e, the average number of hydrogen 
bonds among the water molecules in the monolayer is larger 
than the others, and the number of hydrogen bonds between 
the monolayer and the second layer is the smallest. The total 
number of hydrogen bonds per water molecule in the mono- 
layer is also about 3. Thus, the large charge and the suitable 
unit cell size (J = 0.142nm) make the water molecules in 
the monolayer prefer to form hydrogen bonds within the wa- 
ter monolayer, rather than form hydrogen bonds between the 
monolayer and water molecules in the second layer. Clearly, 
the unit cell size is also the key to the formation of hydrogen 
bonds of the water molecules near the solid surface. 
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IV. CONCLUSION 


In summary, we study the structure, properties of free en- 
ergy and hydrogen bonds of ordered water monolayer on 
ionic model surface with graphene-like hexagonal lattices 
with different charges and unit cell sizes by molecular dy- 
namics simulations. The results indicate that both the charge 
and unit cell size have a great effect on the water molecular 
behaviors in the monolayer, such as water molecular configu- 
rations and the hydrogen bond network. The charged surface 
displaying strong adhesive interaction is described by the wa- 
ter density profiles and potential of mean force. We have also 
carefully investigated the number of hydrogen bonds of the 
ordered water monolayer near the solid surface. It is expected 
that the finding in this paper may help to deeply understand 
the ordered water monolayer on the surface. 
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